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I. INTRODUCTION 


The optically processed Images from a SLAR are useful to Investigators 
In the field of remote sensing to the degree that the wavelength and 
polarization response of terrain and water areas can be quantitatively 
determined. Since the radar video Is recorded on photographic film with 
Its Inherent limitation In dynamic range* It Is desirable to consider 
converting the video to digitally recorded data on magnetic tape* for 
subsequent correlation to r.'^ar maps and analysis of the resolution cell 
contents. The video data collected by the radar are thus preserved In- 
tact and the subsequent operations to recover the Images may be closely 
controlled on a high speed digital computer and high resolution electronic 
display system with controllable dynamic range. The photographic pro- 
cessing* through an analog recording system and an analog optical film 
processor* suffers from the aforementioned limitation In film dynamic 
range (on the order of 23 dB) and correlation errors, which result In 
loss of precision In the photographic Images produced. In addition to 
these errors* which may occur during the correlation process or film 
exposure, tight controls are required on film processing to generate a 
satisfactory photograph. 

With properly digitized data* the precision of the operations can 
be controlled to the level represented by the digital equipment used 1r. 
recovering the data* and the characteristics of the electronic display 
system used. Additionally* the data may be made available on magnetic 
tape records for analysis by users In the field of earth resources. 
Photographs of the display may be used qualitatively to Identify areas 
of interest for detailed evaluation of magnetic tape data. 

The NASA Lyndon B. Johnson Space Center (NASA/ JSC) operates an 
APQ-102 modified SLAR radar on an RB-57 aircraft test bed, normally 
flying at about 60*000 ft* to simulate orbiting vehicle geometries and to 
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gather wide coverage camera data. The work described In this report 
Involved the design, fabrication, and Installation of a video sampling 
and recording system on board the aircraft, and the processor develop- 
ment to generate digital radar lirogery. The work was performed by 
ARL:UT under contracts with NASA/ JSC. 
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II. RADAR AND DIGITAL SAMPLING SYSTEMS 


The modified APQ-102 generates t»M) pulse compressed video signals* 
one from a co>polarized antenna with nominal 3 dB beamwidth of I.S'^* and 
one from a cross-polarized antenna. The radar uses a linear frequency 
modulated pulse with 15 MHz bandwidth at an X-band operating wavelength 
of 3.125 cm. The compressed 3 dB pulse width is about 60 nsec. The 
radar prf is locked to the aircraft ground speed with a sample rate 
of two per foot. The transmitter can be switched to either the 
horizontally or vertically polarized antenna* both of which are gim- 
baled to provide motion compensation and to keep the boresite 
direction normal to the velocity vector at nadir angles of 33 or 54°. 
Figure 1 illustrates the antenna installation. 

The digital data recording system (DDRS) was designed to digitize 
the raw (pulse compressed) video signals from the radar and various 
parameters from the radar operating environment* such as the inertial 
platform signals and aircraft operating parameters. These are 
necessary in reconstituting the radar video Into synthetic aperture 
maps and* subsequently* In analyzing the recorded data. 

The DDRS is installed in the instrumentation pallet which is carried 
In the center section of the RB-57 aircraft test bed. Figure 2 is a 
photograph of the electronics pallet which houses the APQ-102 photographic 
data recorder as well as the DDRS. 

A. Characteristics 

The DDRS receives the analog video signals generated by the 
APQ-102 and normally used to modulate the two axes of two cathode ray 
tubes in the photographic recording system. Buffer amplifiers supply 
both horizontally and vertically polarized, synchronously detected* and 
range pulse compressed signals to the DDRS. The operator can adjust the 
gain of these signals to the DDRS and the photo recorder. In addition* 
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the operator can select the radar mode* which essentially selects the 
dept'esslon angle of the antenna system (or the nadir angle) to change 
the swath of terrain Imaged for a nominal operating altitude of 55,000 
to 60,000 ft MSL. The system also digitizes and records various 
operating parameters of the radar and Its platform, the power supply 
status for various power supply voltages, and operating environment 
temperatures. This permits evaluation of the digitizing and recording 
subsystem In conjunction with each data gathering exercise. 

The details of the 0DR$ subsystem and Its Interface with the 
APQ-102 radar are documented in ARL-TM-80-14,^ an operating and maintenance 
manual produced under NASA Contract NAS9>15217. Figure 3 Illustrates 
the relationships among the various components. 

B. Parameter Options Available 

The principal option available to the system operator Is the 
sampling rate of the DDRS which permits selection of sampling Intervals 
by the high speed analog-to-digital (A/D) converters of 40-150 nsec 
In 10 nsec steps. The system operation at the 40 nsec Interval Is 
critical and requlreii careful setup of the timing pulses for satisfactory 
performance. Coupled with the selection of sampling rate Is the 
precision of the data recorded, I.e., 1-, 2-, 4-, or 7-b1t levels of 
digital data for each sample, which may also be selected by the operator. 
Certain combinations of mode (there are two modes or operating geometries), 
sample Interval (SI), and data precision (DP) are Invalid and are so 
Identified in the recorded data. Finally, the start of sampling Is set 
by the delay selected, 2-99 usec. 

The prf of the radar system Is determined by the ground speed, so 
that two pulses per foot of travel along Its velocity vector are 
generated by the radar. This parameter Is not adjustable, nor Is the 
range resolution, which Is set by the analog pulse compression system 
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of the APQ'102. It has a tine bandnUth product of (1.2 usec)(15 NHz)* 
or a conpresslon ratio of 18:1 1 with a conpressed pulse widtii of about 
60 nsec or 30 ft In slant range. The DORS can thus oversanple at the 
40 and SO nsec SI selections* or undersample at the higher rates witn 
respect to the nominal 3 dB range resolution of the radar. 

The azimuth resolution of a sidelooking SAR Is limited theoretically 
to half the antenna dimension for the focused processing situation; for 
the APQ-102 It Is the antenna azimuth plane dimension of 120 cm/2* or 
60 cm* roughly 2 ft. 

Practically speaking* the resolution In azimutn Is limited to about 
30 ft bv errors Introduced In motion compensation of the video. The 
antenna Is designed to minimize Doppler ambiguities* I.e.* act as a 
spatial filter to eliminate echoes from scatterers with modulo 2« radian 
phase shifts 1<k»nt1ca1 to phase shifts from other (different) scatterers 
In the main beam. This real beamwidth* coupled with the prf* determines 
the niaber of time samples that can be obtained from any given scatterer 
for processing In a syntnetic aperture. The real beamwidth therefore 
sets the maximum length of a synthetic aperture and the resolution or 
main lobe width of the synthetic aperture. 

The main lobe width of the APQ>102 In azimuth Is a nominal 1.3” 

(22.7 mrad) at the 3 d6 points. At the nominal slant ranges for 
operation In modes 1 and 2 at the center of an Imaged swath, the 
alongtrack ground coverages are 1710 and 2478 ft* respectively, so at 
2 samples/ft* a large number ran be processed. To keep the sample 
grid spacing approximately the same 1.: azimuth and range* a much smaller 
number of pulses should be processed In a synthetic array. The choice 
of array length will determine the need for applying a quadratic phase 
correction or focus correction to tl^e recorded data prior to processing 
the synthetic array. 
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An additional option for the processor Is the selection of a 
weighting function, or window, to optimize the tradeoff between width 
of the mam lobe In the Doppler spectrum and the suppression of 
sidelobes. The processor developed by ARL:UT provides either Taylor, 
cosine on a pedestal (Hamming), or Hanning weights for the pulses 
In the arrays. Other weighting routines, such as Kaiser, Chebyshev, 
etc., may be employed by providing the additional algorithms to the 
software. 

Finally, the choice of display parameters must be made In the 
postprocessor, based on the statistics of the Image calculated by the 
program from the overlaid outputs of the PoppIPf* filters or pixel 
values. The operator must also choose the threshold of the displayed 
data, dynamic range, and logarithmic or linear data options, as well 
as the size of the displayed coverage. The latter is fixed In the 
ARL:UT processor at 360 azimuth lines by 384 range sample intervals. 

In essence, there'iore, one set of parameters is chocen at the 
equipment, and another set in processing the raw data. Table 1 provides 
a breakdown of these choices. 


TABLE ! 

SYSTEM SELECTABLE PARAMETERS 


Operator Selectable Parameters 

DDRS Radar 


DRMIN - Delay time to start of sampling 
SI - Sampling interval for A/D 
DP - Data Precision (1 ,2,4,7 bit) 


Mode (1 or 2) 

Gain (H or V video) 
Transmitted Polarization 


Note: The DDRS design incorporates logic to flag incompatible 

or invalid combinations of the above parameters. Valid settings 
are listed in Table C-1 of ARL-TM-80-14 for various combinations of 
mode, delay, SI, and DP. Nevertheless, data can still be re- 
corded, despite invalid control settings. 


Processor Selectable Parameters 
Co-polarized or Cross-polarized Data 

Area to be Imaged: Starting Range Sample + 384 Range Bins 

Starting Frame Number + 360 Azimuth Elements 

Azimuth Resolution - Fixed by the Software to be the Same as 

Range Resolution at the Map Center 

Array Weighting and Beam Broadening Factor 

Range Bin Used as the Map Center 

Range Bin Used as the First Range Bin on the Image 

Postprocessor/Display Parameters 
Log or Linear Data 

dB/Gray Shade or Number of Filter Magnitudes/Gray Shade Assignment 
Brightness, Contrast Settings 
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III. RECORDED DATA 


The DDRS was designed to record the digital data, video, and radar 
parameters on a 14-track wideband tape recorder in Miller code. The 
format was designed to be compatible with the DRIE equipment used at 
ARLrUT. 

A. Formatting and Track Assignments 

Details of the data formats and the track assignments on the Mincom 
(in-flight) tape recorder may be found in ARL-TM-80-14.^ Unfortunately, 
tho track assignment utilized by the NASA/JSC contractor does not agree 
with that prescribed, requiring some cabling adjustments in the ARL:UT 
playback equipment for compatibility with the reformatting hardware. 

B. Decoding and Reformatting Requirements 

The video data and the radar parameter environment (RPE) data are 
recorded on the wideband tape after being Miller encoded. The Miller 
code has certain advantages over the NRZ (non-return-to-zero) code in 
recording efficiency and S/N ratio. Miller decoders are required for the 
8 tracks of video (4 horizontal and 4 vertical polarization) and 1 track 
of RPE data. So-called NERDAS (NASA Earth Resources Data Acquisition 
System) data and time codes are recorded directly in NRZ code. The 
digital video and RPE data, after decoding, are processed through the 
DRIE to an IBM computer-compatible format and rerecorded on a Hewlett- 
Packard 7970B tape drive unit. Figure 4 illustrates the process. The 
DfilE reformats the data as described in the technical report ARL-TR-76-42. 
The wideband tape formats are described in the technical manual 
ARL-TM-80-14.^ 
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C. DRIE - Functions and Characteristics 


The ORIE equipment was designed for the Flf,MR radar system, and 
hence Its application to the APQ-102/DDRS Is not very efficient. It 
duplicates a realtime airborne processor, and the function of converting 
the data. Illustrated in Fig. 4, Is not Implemented adequately for this 
purpose. Depending on the bit precision of the data recorded, a maximum 
of 128 range bln samples may be reduced at one pass through the selected 
frame numbers (FN) on the wideband tape. This Is true for the 4-, 2-, 
and 1-b1t data. With 7-b1t data, however, only 16 range samples can be 
processed at a time; consequently, 24 passes are required through the 
wideband tape records to convert 384 range samples of data. Usually, 8-12 
9- track, 800 bpi computer tapes are required for a full image set of raw 
data. 

The ORIE contains bit parity error detecting circuitry which permits 
monitoring the source and rate of bit error generation, a feature of 
considerable value in evaluating the performance of the DDRS and wideband 
tape recording equipment. It also indicates error sources in the re-encoding 
process, such as the Hewlett-Packard tape drive buffer memory and tape 
recorder operation. Various addresses for the RPE data may also be dis- 
played, with their bit contents, to evaluate the data as decoded from the 
instrumentation (wideband) tape. 
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IV. PROCESSING REQUIREMENTS 


The video data are recorded in continuous strips. The range swath 
beginning and width are determined by the ORMIN variable (the selected 
delay following the STC trigger in the APQ-102), the sampling interval, 
bit precision, and mode (see Performance Specifications, pp. 20-21, 
ARL-TM-80-14^). The range coverage is thus determined by the number of 
samples to be displayed and the sample interval. Likewise, the displayed 
azimuth coverage is determined by the Doppler filter separation or 
spatial frequency corresponding to spectral frequency, and the number 
of unique filters, after overlay. The ground signal processor (GSP), 
in its delivered configuration, makes the azimuth interval equal the 
range interval at the center of the image. 

A. Physical Limitation on Parameter Selection 

The GSP options are limited by the previously selected values of 
the DDRS and APQ-102 parameters, and by certain design parameters of the 
radar system. The most important of the latter is the fixed relationship 
between the platform ground speed and the radar prf of 2 pulses/ft, 
which sets the maximum synthetic beamwidth of the system and the fixed 
slant range resolution capability of the radar (about 30 ft). Sampling 
in range can be varied from a 1.5 oversampling ratio at the 40 nsec rate 
to 0.4 undersampling at the 150 nsec rate. Once the DDRS parameters are 
set by the airborne operator, the GSP characteristics are fixed, except 
those that determine the azimuth resolution and at the same time, the 
sampling ratio, since the sampling rate is fixed by the platform velocity. 
Second, the GSP operator must select the portion of the sampled radar 
data to be processed, depending on the number of display elements available. 



The smallest number of range samples recorded (at 7-bit precision) 
is 512, and this is larger than, for example, the number of display 
elements (pixels) on the ARL:UT High Resolution Display System. The 
number of azimuth elements is usually of the same order as the range ele- 
ments on most CRT displays. Thus, the displayed area should be small 
enough to accept (digitally) a legend containing a description of the 
scene, its dimensions, and data concerning the flight parameters. 

Finally, the GSP operator must select the postprocessing display 
parameters such as the dynamic range represented, number of bits, 
intensity resolution, etc. This selection can be made after examining 
the statistics of the linear data for the radar image, i.e., the mean, 
standard deviation, and maximum and minimum values generated after 
array processing and overlaying. The characteristics of a display will 
determine this selection, i.e., number of gray shades, availability of 
color, etc. The ARL:UT display is limited to 16 gray shades, and the 
dynamic range is determined by the assignment of dB/gray shade. Usually, 
2 dB/gray shade is used, with a dynamic range of 32 dB. 

B. Processor Organization 

Figure 5 is a flow diagram of the processor, which is organized into 
a Doppler processor and a postprocessor. Table II lists the routines 
included in the GSP. 

The first step is to reformat the DRIE data records. Program 
REFORMP reads the DRIE records, decodes them, and then outputs them to 
a scratch disk file. REFORMP also has the option of checking the DRIE 
records for errors; each DRIE record contains parameters which alert 
the operator to possible errors in the radar, DDRS, and DRIE hardware 
and its tape output unit. The first file generated by REFORMP contains 
recorded radar parameters. The rest of the generated files contain 
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TABLE II 

PROGRAMS AND SUBROUTINES INCLUDED IN THE GROUND SIGNAL PROCESSOR 


GSP 

REFORM? 

REFORMD 

RDECOD 

RSCALE 

HVDECOD 

IBCDTD 

APWT 

TAYLOR 

PULSES 

LIGHTF 

LUNPOS 

FFTCC 

MAPSTAT 

DSTATR 

PIMAGE 

ORTA 


Main program, performs Doppler processing and filter 
overlay 

Program, decodes DRIE records and formats them into 
range bin records 

Program, uses output of REFORMP and creates files of 
samples by range bin 

Subroutine used by REFORMP to decode the RIOT parameters 
from DRIE formatted records 

Subroutine used by REFORMP to interpret the decoded RIOT 
parameters from part of the output of RDECOD 

Subroutine which recovers the video data samples from 
DRIE formatted records 

Function, converts BCD to binary format, used by RSCALE 

Subroutine used to calculate an array of weights for 
different aperture weighting functions 

Subroutine used by APWT to calculate an array of weighting 
coefficients for the Taylor weighting function 

Subroutine, it skips and recovers video pulses from the 
input data file into an array for use by GSP 

Function which checks the status of pseudo-sense switches 
set by the operator before program execution 

Subroutine that positions data files 

Subroutine, computes fast Fourier transform of a complex 
sequence; from the International Mathematical and 
Statistical Library 

Program which reads all or part of GSP output files of 
filter magnitude values and compiles statistics on the data 

Subroutine which compiles and computes statistics, used by 
MAPSTAT 

Program, performs the postprocessing on filter magnitude 
data in files generated by GSP, then outputs the data in 
records formatted for input to the ARL:UT high resolution 
display software 

Subroutine, reformats azimuth lines of 384 range samples 
into records compatible with the ARL:UT high resolution 
display software, used by PIMAGE 
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1024-word records of sampled video data at different range bins. Program 
REFORMO Is run next to reorganize the records into files containing the 
sampled video at Individual range bins. At the start of the reformatting 
process, the operator must choose to decode either the co-polarized or 
the cross-polarized data. 

The output file from REFORMD Is used as Input to the main program, 

6SP. GSP requires the operator to select parameters such as the aperture 
weighting function, shading factor and, *n particular, the range bln of 
the map center, an Important parameter In the processor algorithms. 

Certain default values are provided In the programs. For example, 
the Taylor aperture weighting function with 30 dB peak si delobes Is the 
default for the window function, and 1.27 (BBF) is the 3 dB width of the 
resulting Doppler filters. If a different weighting function is chosen, 
the BBF must be calculated or obtained, and entered into the program 
as a change. 

The selection of synthetic array length is governed by the range 
sampling Interval originally selected and the operating mode and patch 
center, which determine the nadir angle to the patch center. Referring 
to the GSP program listing, the Doppler spectrum line spacing Is set 
equal to the projection of the slant range sample spacing on the horizontal 
or Image plane. This Is fixed In the software, but could be changed by 
making FILTSPC equal any arbitrary spacing. 

In using digital SAR data to provide radar Images for remote 
sensing applications, achieving maximum resolution In azimuth was not 
considered Important. In addition, by keeping the array length shorter, 
focusing requirements are less stringent, although focusing (quadratic 
phase) corrections are Included. 
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Following establishment of the window weights, the focus weights 
are obtained and placed In a complex array FOC. The correction to the 
1th pulse In the array Is 

OPHI » - * -~(VRANGE^ - range) , 

where 0 is the distance from the 1th pulse to the center pulse In the 
synthetic array. The phase angles are converted to 1n»phase and 
quadrature components and combined with the aperture weights. 

Since the ground coverage of the synthetic arrays Is a function 
of range, the corresponding filters from arrays at different range bins 
may not i^epresent the same azimuth components. To maintain the integrity 
of the image, the difference In the azimuth coverage at the current range 
bin versus the coverage at an initial reference range bin Is computed. 
This reference range bin is input as parameter IRBFMAP. If processing 
Is to begin with range bin one, then IRBFMAP*!. The azimuth component 
represented by the first filte** from the first available synthetic array 
of range bin IRBFMAP is then the azimuth coordinate of the start of the 
map. To have the first filter of following range bins represent the same 
azimuth coordinate, the difference in the coverage Is converted into 
pulses and these pulses (data samples) are skipped before Doppler 
processing begins. 

Subroutine PULSES skips data samples In a file and recovers data 
Into array VIDATA. The data for one array are then weighted with the 
proper weight from array FOC and placed in a complex array HV. The 
array is processed by a complex FFT and the resultant filters are 
obtained by summing the square of the real and imaginary components 
out of the FFT. The corresponding filters from consecutive arrays are 
summed (overlaid), and the number of overlays at each azimuth coordinate 
is kept In array RNORM so that the filters may be normalized for the 
number of summations they contain. 
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C. Image Formation 


Referring to Attachment I for the algorithms. Attachment II for 
numerical values associated with the Images shown in Section V, 

Table III for a list of 6SP variables, and Fig. 6 for the flowchart, 

rESR . CL * SISR . (0.2998x10^1(1. 006x10'^ 

2 * SIN(XNAPC) 2 s1n(0. 70973 

RESR » 23.14 m and FILTSPC is set to this value. 

Each input file into 6SP contains the samples from one range bin. 
The coverage and filter spacing of the synthetic arrays are a function 
of range. Based on a desired filter spacing on the ground of FILTSPC, 
the size of the synthetic arrays is calculated, 

NP (at range R) » p^LTsPC * ' ‘ 

Since NP is an integer, the filter spacing out of the synthetic arrays 
is actually given by FILTSP, 

FILTSP » • 

This number is calculated and printed for reference; however, for the 
parameters encountered in this system, FILTSP and FILTSPC are nearly 
equal and FILTSPC is assumed to be the filter spacing. 

The 3 dB resolution of the synthetic filters is also calculated 
and given for reference, 

«... _ 0.4429 * BBF * WL * RANGE 
“5A (NP-1 ) • OTH ■ 

Next, the pulse weights for the input aperture weighting function are 
calculated in APWT. 
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TABLE III 

VARIABLES AND CONSTANTS USED BY PROGRAM GSP 


PI 

TP I 

PIH 

RTD 

DTR 

FPM 

CL 

WL 

BBF 

KWTFTNA 

SHADFAC 

NRBPC 

IRBFMAP 

DRMIN 

SI5R 

NBIT 

NRB 

RMIN 

RPC 

XNAPC 

IRBF 


The value of ’f 

Twice the value of v 

Half the value of n 

Degrees/radian, conversion factor 

Radians/degree, conversion factor 

Feet/meter, conversion factor 

Velocity of light, 2.997925E3 

Operating wavelength of radar, 0.031228 m 

Beam broadening factor, default»1.27 

Name of weighting function in H format, default*''TAYLOR" 

Scale factor for the weighting function; for Taylor 
weighting, it is the peak sidelobe ratio (In dB), default>30 

Number of range bins to be used as the map center, 
default“l/2 of the number of range samples per pulse 
taken by the DDRS 

Number of range bins to be used to set the azimuth 
coordinate of the start of the image, default*range bin 
of the first file processed 

Delay (sec) from the S7C trigger to the start of 
sampling by the DDRS 

Slant range sample interval of DDRS 

Bit precision of the sampled video data 

Number of digitized samples or range bins, per pulse 

Slant range distance (ra) from the radar to start of 
DDRS sampling 

Range to the map center, m 

Nadir angle to the map center 

Range bin of first data file processed 
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TABLE III (Conf d) 

VARIABLES AND CONSTANTS USED BY PROGRAM GSP 


RNGFMAP 

RESR 

FILTSPC 

NPT 

lAZLINE 

NRBPROC 

RANGE 

XNP 

NP 

NF 

FILTSP 

RESA 

WT 

GI 

DPHI 

FOC 

NPLMAP 


VI DATA 

HV 

IWK 


Slant range distance, m, from radar to the range 
bln IRBFMAP 

Ground spacing, m, between range samples at the map 
center 

Desired spacing (ft) of output azimuth filters 

Counter to keep track of total number of pulses used 
during processing of each range bin file 

Counter pointing to the azimuth line with which the 
first filter out of the current array is overlaid 

Range bin currently being processed 

Slant range (m) to the range bin currently being 
processed 

Number of pulses to process per array, real variable 
Round of XNP 

Number of filters recovered per synthetic array 

Actual filter spacing (ft) of the current arrays 

The 3 dB resolution of the Doppler filters of the 
current arrays 

Arrry containing the amplitude weight for each pulse 
in a synthetic array 

Inverse gain factor of the weighting array, from APWT 

Quadratic phase correction angle 

Array containing the focusing corrections and 
amplitude weights for each pulse 

Number of pulses to skip on the current range bin file 
so the first filter out of the first synthetic array 
represents the azimuth coordinate of the left edge of 
the map 

Array containing the video data samples of one synthetic array 
Complex array used for input and output to the FFT routine 
Integer scratch array for FFTCC 
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TABLE III (Cont'd) 

VARIABLES AND CONSTANTS USED BY PROGRAM GSP 


WK 

RMAP 

RNORM 

DFLOMN 

NPSKIP 


Scratch array for FFTCC 

Array containing the processed overlaid filter 
magnitude values for a range bin 

Array containing the number of overlays at each 
azimuth line 

Distance (ft) the radar has moved since the start 
of azimuth processing 

Number of pulses to skip before forming the next 
synthetic array 
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Based on the principle of keeping the Doppler filters aligned 
from one array to the next* It is necessary to discard some of the 
pulse samples. The distance flown (In feet) Is computed as equal to 
half the total number of pulses* NPT. reached at the last pulse of 
the last array processed, based on two pulses/ft. Next, the residual 
fraction of a filter spacing In azimuth Is calculated by dividing 
OFLOWN by FILTSPC and discarding the Integer part of the quotient. 

The fractional part. FLEFT, Is subtracted from a filter spacing and 
the remainder Is converted to an Integer number of puUes to be skipped 
by doubling, adding 0.5, and truncating. 

NPSKIP = 2(RESAEF-FLEFT) + 0.5 

The next pulse (after skipping) becomes the first in the next array. 
Finally, the overlaid filter magnitude values in RMAP are normalized, 
using RNORM, the number of overlays. These data are buffered out onto 
tape or disc file. 

The user must control the processing by interacting with the 
processor to provide various parameter or data selections. 

The principal operations consist of the following. 

1. Select the desired frame numbers (octal) to start and stop 
the map. The start FN is usually the first one on the file. 

2. Indicate the tape number and select the file to be processed, 
which usually contains a strip of 16 or 32 range samples. 

3. Select the polarization you wish to process. 

4. Select the aperture weighting function and shading factor. 
Default Is Taylor weighting, 30 dB PSLR. 

5. Select patch (map) center range bln. 

6. Select beam broadening factor. For Taylor 30 dB weighting, it 

2 

is 1.27. Cosine and Hanning must be determined after establishing 
shading factor. 
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The main processor program GSP, listed in Table II» is supported 
by a number of subroutines and library functions. The most important 
ones are HVOECOO* which decodes the horizontal or vertically polarized 
video signals; RDECOD, which decodes and makes available necessary 
radar system operating parameters; and a complex Fourier transform 
routine FFTCC from International Mathematical and Statistical Library. 
Others include APUT, the calculation of synthetic aperture pulse 
weighting values; PULSES and RECOVRY. used in organizing and obtaining 
the data from tape and disc; and LIGHTF, a means for checking a pseudo- 
sense switch. Program MAPSTAT compiles statistics, i.e., minimum and 
maximum values, mean and standard deviation. MAPSTAT may be run on the 
entire map or any desired section. Figure 6 is a flow diagram of GSP. 

D. Synthetic Array Coverage 

Consider a digital time series of samples of video data in a 
synthetic array. The sample rate is the PRF of the radar so that PRF/2 
is the highest recoverable frequency in the series and Af=PRF/(Np-l ) is 
the spacing between resultant filters of a Fourier transform performed 
on the series. 

The two-way Doppler shift associated with a target is given by 

.. . . 2V sine 

■ "x 1 • 

where is the radial velocity, X is the operating wavelength of the 
radar, and e is the angle from the normal to the target. Figure 7 
illustrates the geometry of the ground coverage of a synthetic array; 
for typical values of R, i.e., slant range from radar to ground, a«3° 
so that sinaiwa«x/R. As stated above, the maximum recoverable frequency 
is PRF/2 so that 

^f = PRL « 21 sing , m 
d 2 X ^ 
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FIGURE 7 

SYNTHETIC ARRAY COVERAGE GEOMETRY 


ARL;UT 
AE*81-87 
CRG • GA 
5-21 -81 




Nom, two pulses per foot are emitted by the radar, so PRF*2V If PRF 
Is In hertz and V In ft/$ec, and x»Ra/ 2 If R, a are In feet. 

Thus the coverage for any synthetic array formed Is a function of 
the range represented by the data samples of the array. The radar Is 
mounted on the left side of the plane so the Image Is made from left to 
right as the synthetic arrays are formed. To start the coverage of 
the arrays at the same azimuth coordinate, the difference in coverage, 
A, at R and R^ (RNGFMAP In program 6SP) Is computed. Then a ft x 2 
pulses/ft Is the number of pulses to skip before processing the range 
bln at R. 


The spacing between the resultant FFT filters, d^, represented on 
the ground is 

^f " (Np^T . R In ft . 

Thus, once a value for the filter spacing Is obtained (In GSP It Is 
arbitrarily taken to be the spacing between range samples at the patch 
center), the number of pulses, Np, necessary to obtain that spacing Is 
fixed. In this way the spacing between azimuth filters remains the 
same at different ranges. Since Np must be an Integer, one can only 
get close to the desired spacing; however, with the parameters and 
geometry encountered In this system, actual spacings within 1% of the 
desired spacing are realized. 


Finally, the 3 dB resolution of the formed synthetic beams (filters) 
Is given by 


0.4429AK^R 

Pa “ (N -1)0.5 ft/pulse 
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where K Is the beam broadening factor resulting from nonuniform 
weighting of the synthetic arrays. Note that “ 0.886 K^d^; thus 
for Taylor weighting (K_*1.27) the 3 dB resolution of the synthetic 

fl 

beams Is approximately 1.12 d^. The only variables which affect the 
accuracy of the calculation for filter spacing are x. R. and the 
accuracy of PRF versus ground speed, I.e., how accurate Is the 
assumption of two pulses per foot emitted by the radar. Inaccurate 
determination of d^ results In Improper overlay, which has the effect 
of degrading the azimuth resolution of the Image. 


V. POSTPROCESSING AND DISPLAY 


Attachment III is a listing of the short program which generates 
the image display tape. In generating the tape, the operator should 
refer to the (optional) printout of statistics of the filter magnitude 
values. This will provide insight to the gray shade resolution, noise 
threshold, and dynamic range of data to be displayed. 

A. Display Parameter Options 

There are three display parameter options: 

(1) Logarithmic or linear data. Logarithmic data are most 
representative since the human eye functions in a logarithmic manner, 
and dynamic range is compressed on the display. 

(2) Gray shade level assignments. The gray shade resolution may 
be adjusted by selection of the levels per gray shade, DB6S assignment. 

A default value of 2 is provided in the software. 

Note: Dynamic range is set by (2) above. The ARL:UT high 

resolution display has 16 gray shade levels, so DBGS x 15 is the dynamic 
range . 

(3) Threshold. A noise cutoff threshold is set by the assignment 
of a selected level to the top gray shade (number 15 out of 0-15), and 
then subtracting the dynamic range. 

Once the data are displayed, the characteristics of the display 
itself must be considered. In particular, the brightness and contrast 
controls will affect the quality of the displayed radar image. 
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B. Image Quality Considerations 


In providing maximum Image quality, a most Important consideration 
Is resolution. In this case, the sampling rate of the DORS sets the 
resolution of the Image. Comparisons of digitally and photographically 
processed Images reveal the effects of undersampling. Figures 8 and 9 
illustrate the evolution of high quality digital Imagery. 

1. Radar/DDRS Parameter Selection n^^uons - Selection of the 
longer SI (sampling Intervals) on the DORS increases the range swath 
recorded, but results in undersampling in both range and azimuth 
(although the latter could be changed in the GSP). The choice of bit 
precision also affects the quality of the image, since the lower bit 
precision values result In decreased data resolution compared with 
7-bit data. Mode selection affects ground resolution versus slant 
range resolution, since In mode 2 the cosecant of the nadir angle is 
less than in mode 1. 

2. Processing Options - Since the factors above are set, once 
the data are recorded, only the selection of polarization channel to 
be processed and those factors discussed previously (such as the 
weighting function), can be adjusted. In particular, the display 
options in Section V.A become effective. Figures 8 and 9 illustrate, 
to some degree, variations in gray shade assignment, focusing of the 
data in a relatively short array, and overlaying of arrays. These 
variations affect the displayed image quality in varying degrees. 

Although the processor does not generate images with maximum azimuth 
resolution, the saving in processing time was considered the overriding 
factor. If maximum length arrays were processed, then the focus function 
would become very important, whereas Improvement over the unfocused 
image for the relatively short arrays is almost indiscernible. 
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C. Algorithm Details 


Procram PIMAGE is the postprocessing program which requires the 
inputs listed in Section V.A. After entry of DB6S and T15, (power 
level assigned to the top gray shade) a threshold is set, THOLD » 10 
log(T15)-16 DBGS. Then an integer value from 0-15(IGS) is assigned 
to each p.xel power level (Pi) output from GSP, 

IGS * (10 log(Pi)-THOLD)/DBGS 
If IGS<0, then IGS=0; or if IGS>15, then IGS*15. 

Next, subroutine ORTA is called to output the display values 
IGS in a format compatible with the ARL:UT high resolution display 
software. Following tne same procedure, values could be assigned 
based cn any number of gray shades and output in formats compatible 
with other displays. 
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ATTACHMENT I 
LISTING OF PROGRAM GSP 
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ATTACHMENT III 
LISTING OF PROGRAM PIMAGE 
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